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Recently, Ding et a/. 1 investigated spin-polarized 
transport in graphene irradiated by a linearly polarized 
laser field. There are several serious problems in their 
model, such as the violation of the charge conservation 
in the graphene region, the incorrect application of the 
rotating- wave approximation (RWA) and even the wrong 
Green's functions. In the following, we discuss these 
problems point by point. 

We first demonstrate that their approach violates the 
charge conservation in the graphene region. Our starting 
point is Eq. (24) in Ref. [J which gives the time-averaged 
current flowing into the left lead. That equation reads 



(II) = 
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where r =t (I) represents the spin-up (-down) band; 
a = L (R) stands for the left (right) lead; the chemical 
potential in the leads are (j*l,r — ±eV/2. For conve- 
nience, we transform all relevant Green's functions and 
self-energies from the rotating reference frame into the 
fixed reference frame, i.e., 
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with s — 1 (—1) standing for the conduction (valence) 
band. Equation (fT]) is then rewritten as 
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in which 
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with = rj^ = r (l ±p). The Green's functions 
glZ' r (s) are given as Eqs. (A.3)-(A.8) by Ding et al. in 



Ref. [l|. However, their expressions are mathematically 
incorrect, particularly the terms given by Eq. (A. 8) are 
even dimensionally incorrect. We present the correct for- 
mulae of gls'' T ( £ ) m Appendix IA1 in order to avoid mis- 
leading. 

Exchanging L and R in Eq. ([3]) , one obtains the current 
flowing into the right lead, 
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Thus 
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with G^ m {s) = EC, SS '(4 To show F T (e) analyti- 
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cally, we first neglect the interband term (i.e., s ^ s') of 



0.06 



0.04 - 



0.02 ■ 



with interband terms 
without interband terms 
field-free 



-2 

e(o) ) 



FIG. 1: (Color online) F r defined by Eq. © as function 
of energy e with and without the interband term of g^[(e) 
in the case with the field strength Eo = 1200 kV/cm and 
frequency ujq = 0.04t 9 . The other parameters are p — 0.4, 
To = 0.05t g , V g = and D = 3t g , same as those used in 
Fig. 1 in Ref. [J We also plot the results in the field-free case 
as the green dotted curve. 
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with D T (e) = — -^Im 9 oo. ss( £ ) being the density of 

states of the isolated graphene under the laser field. 
Since D T (e) is positive semidefinite, F T {e) is also posi- 
tive semidefinite. Also considering / Q (e) > 0, one obtains 
dNo/dt > 0, indicating the accumulation of the charge 
in the graphene region. 

To show the influence of the interband term, we plot 
the numerical results of F*(e) with and without the in- 
terband term of g^mi 8 ) m Fig- HP Here we use the field 
strength E n = 1200 kV/cm and frequency lj = 0.04i s . 3 
The other parameters are p = 0.4, T = 0.05t g , V g = 
and D — 3t g , same as those used in Fig. 1 in Ref. [l|. 
From this figure, one observes that F^(e) with and with- 
out the interband term of g^mi 5 ) almost coincide and are 
both positive semidefinite. One clearly concludes that 
the approach reported by Ding et al. violates the charge 
conservation in the graphene region. Even worse, this 
problem not only appears in the laser-applied case dis- 
cussed in Ref. [l|, but also appears in the field-free case^ 
[^(e) in that case are plotted as green dotted curve in 
Fig. HJa)] and the case with the time-alternating gate 
voltage. 6 Therefore, all results in these works are scien- 
tifically incorrect Z 

Another severe problem in Ref. [l] is that they mis- 
takenly applied the RWA to the case with strong laser 
field in the whole momentum regime. As shown in our 
recent work^ the RWA is only valid for the weak laser 
field at the momentum around the resonant point, i.e., 
2vpk = loq with ujq being the frequency of the laser 
field. In order to make this issue more pronounced, we 
plot the sideband quasi-energies and weights [defined by 
Eqs. (10) and (11) in Ref.||] with and without the RWA in 
Figs.[2ja) and (b), respectively. The parameters are the 
same as those of Fig. [TJ Here we only show the case with 
the momentum along the current direction, i.e., the di- 
rection along the x axis, as Ding et al.X The Hamiltonian 
with and without the RWA are given by Eqs. (10) and (7), 
respectively in Ref. [J The corresponding eigenstates in 
these two cases can be obtained via the standard Floquet- 
Fourier approach widely used in the literature^— Com- 
paring Figs.rSJa) and (b), one finds that the quasi-energy 
spectrum under the RWA is qualitatively different from 
the exact one. In particular, a huge gap opens around 
the Dirac point in the quasi-energy spectrum under the 
RWA, in consistence with the gap in the bias dependence 
of the differential conductance in Ref. [l|. However, this 
gap is absent in the exact quasi-energy spectrum, as re- 
ported in the previous investigations on graphene under 
a linearly polarized laseri 8 i 12 i 13 

Although the pronounced discrepancy in the quasi- 
energy spectrum with and without the RWA is a con- 
vincing evidence of the invalidity of the RWA in their 
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FIG. 2: (Color online) Quasi-energies of the sidebands £JL 
against the normalized momentum with (a) and without (b) 
the RWA. The color coding represents the weight WJL of the 
corresponding sideband (note that it is a dimensionless quan- 
tity). The parameters are the same as those of Fig. [T] 



cases, in order to nail down this issue, we further demon- 
strate that even under their framework, the differen- 
tial conductance with and without the RWA are qual- 
itatively different. Our approach is as follows. By ex- 
ploiting the eigenstates obtained above, which have the 



form |$ k ^W) 



E 



\4>l 



one obtains 



the Green's function of graphene without connecting the 
leads M 
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We then calculate the Green's function of graphene con- 
nected with the leads via Eq. (U). One obtains the cur- 
rent via Eq. ([3]) and the differential conductance G — 
—d{Ii) /dV . We have examined that this approach can 
recover the results under the RWA from the approach 
of Ding et al. after correcting the errors in their Green's 
functions [Eqs. (A.3)-(A.8)] in that paper (see the correct 
formulae in Appendix |A")) . 15 The conductances with and 
without the RWA are plotted in Fig. [3] It is seen very 
clearly that the conductances in these two cases are qual- 
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FIG. 3: (Color online) The differential conductance as func- 
tion of bias with and without the RWA. The parameters are 
the same as those of Fig. [1] We also plot the results in the 
field-free case as the green dotted curve. 



itatively different, especially the pronounced gap around 
the Dirac point does not appear in the conductance with- 
out the RWA, in consistence with the behaviour of the 
quasi-energy spectrum. Beyond all doubt, the above re- 
sults can be seen as a smoking gun of the invalidity of the 
RWA in the cases discussed in Ref. [l]. Therefore, their 
main results, especially the pronounced conductance gap 
around the Dirac point, are scientifically incorrect. 

This work was supported by the National Ba- 
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No. XDB01000000. 



Appendix A: Green's functions under RWA 

In Ref. Ill, the Green's functions of graphene without 
connecting the leads g T Jf{ £ ) given by Eqs. (A.3)-(A.8) 
are not the correct solution of Eq. (A. 2) in that paper. In 
the following, we present the correct formulae of g T s r Ji (s) 
in the cases with V g — 0. 
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with S representing the area of sample. v v = ^i 2 £ ■ 
vsj\e 2 - A 2 | with v = 0, ± and 
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Comparing the above formulae with Eqs. (A.3)-(A.8) 
in Ref. [l|, one observes a lot of differences between them. 
Specifically, in g[l' r and gJJ""' they missed ujq/2 in all v v 
terms. The errors in the interband term are even more 
serious: their formula can be written as 
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It is seen that all terms in the above equation for |e| > A 
take the incorrect dimension. Since there are so many 
errors in their formulae, these errors are not likely to 
only come from typos. 



4 



In order to show their errors more clearly, we calculate 
the DOS of the isolated graphene in the field-free case, 
which has the form 



DS(s) = -ilmE« ss (£ 

s 

1 ^ Tr ( UJq \ 



A=0 



A=0 



From our formulae, one obtains 
S 
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s\e(\e\<D), 



(A8) 
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which is exactly the well-known formula of the DOS in 
graphene^ and is linear with energy. In contrast, from 
the formulae by Ding et ai, one obtains 
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Obviously, the above formula is incorrect, especially it 
gives the negative DOS when |e| < 
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